arXiv:1505.01595v2 [cond-mat.mes-hall] 29 Jan 2016 


Bloch states in light transport through a perforated 

metal 

Zhyrair Gevorkian a,b , Vladimir Gasparian 
California State University, Bakersfield, California 93311-1022, USA 
Emilio Cuevas 

Departamento de Fisica, Universidad de Murcia, E-30071 Murcia, Spain 

a Yerevan Physics Institute,Alikhanian Brothers St. 2,0036 Yerevan, Armenia. 
b Institute of Radiophysics and Electronics,Ashtarak-2,0203, Armenia. 


Abstract 

Light transport in a metal with two-dimensional hole arrays is considered. An¬ 
alytical expression for a transmission coefficient in periodic, isolated and disor¬ 
dered cases are derived, assuming the existence of transverse waveguide modes 
tunneling in two-dimensional plane perpendicular to traveling direction of light. 
The one dimensional case of periodic holes, due to its simplicity, is investigated 
in detail. In the dilute metal regime, when metal fraction is small, our numerical 
study of the transmission coefficient of central diffracted wave indicates the exis¬ 
tence of a minimum which is completely independent of an incident wavelength. 
Further increasing of metal fraction leads to the unusual monotonic increasing 
of central diffracted wave transmission. The role of the surface plasmons is 
discussed. 

Keywords: periodical hole arrays, Bloch states 


1. Introduction 

Since its discovery [l] the extraordinary optical transmittance (EOT) at¬ 
tracts great interest. This interest is largely motivated by a recent progress in 
nanotechnology which allows to get for EOT a possible applications in different 
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optical devices. Many experimental and theoretical papers devoted to a study of 
the EOT (for a review, see Ref. Q). The phenomenon of EOT seemed to be well 


understood (3, [4j with the involving of the surface plasmons and Bloch states. 
The former appeared on the interface between metal and dielectric, while the 
Bloch states originating by the periodicity of hole arrays. However, so far in 
theoretical understanding of EOT much attention was paid to plasmon aspect 


of the problem 


and less attention to the Bloch states and periodicity. 


The Bloch states of plasmons on the periodically perforated metal surface were 


studied in 7]. It has shown, that EOT phenomenon can be explained using the 
mechanism of plasmons’ vertical tunneling from one metal surface to another 
one, where plasmons are eventually converted to free photon states. 

In the present paper, without pretending to a complete mathematical de¬ 
scription the theory of EOT, we develop a different approach to study the 
behavior of electromagnetic waves in a periodically performed metal system 
taking into account the existence of transverse waveguide modes tunneling in 
two-dimensional surface perpendicular to traveling direction of light. We have 
concentrated our attention on the role of Bloch states and show that this mech¬ 
anism with a transverse waveguide modes tunneling leads to peculiarities in 
light transport. Particularly, in the dilute case when metal fraction is small, 
we have found an independence of transmission coefficient of central diffracted 
wave on the incident wavelength and its unusual increase with increasing of a 
metal fraction. We explicitly derived an analytical expression for transmission 
coefficient in disordered hole arrays case. We found that without Bloch states 
the disordered hole arrays lead to broadening of spectral shape in accordance 
with the Refs. 


2. Formulation of the problem 

Let us consider a metallic film with periodic array of two dimensional holes 
(see Fig.l). Suppose a plane wave enters the metal film from the 2 < 0 half¬ 
space at normal incidence. In order to find the transmission amplitude , we 
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Figure 1: Geometry of the problem. 


start from a scalar Helmholtz wave equation 

V 2 $(x, y, z) + k$e(x, y)<Z>{x, y , z ) = 0, (1) 


where fco = lo/c is the wave number corresponding to the angular frequency 
u> of an incident photon and e( x, y) is the two dimensional periodic dielectric 
permittivity of the system. Eq. m is valid for s and p-polarized waves and 
the scalar function >!> describes the transverse components of an electric or a 
magnetic fields, respectively. 

Mainly following Refs. 10, [ll|, we seek the solution of the propagating in 
the system wave as a product of a fast and a slowly varying, y , z), function 
on a wave incident z direction, implying that the metal fraction of our system 
is small. 

$(x,y,z) = e lkoZ (j)(x,y,z). (2) 


Note, that this situation differs from the most of EOT considerations where 
usually the opposite case is considered. 

Substituting Eq. © into Eq. © and neglecting the second derivative of (j> 
with respect to z (\d 2 (j)/dz 2 \ « 2ko\d(/)/dz\), one gets 

= H{x,y)(j>, (3) 

where 

H(x,y) = ~^ V ? + Y (1-^2/)) ( 4 ) 

and Vj = (d 2 /dx 2 + d 2 /dy 2 ). 

The obvious similarity of Eq. (0 (the spatial coordinate z plays the role of 
the time) and the time-dependent Scrodinger equation for a particle with mass 
fco, moving in the two-dimensional potential V(x,y) = (1 — e(x, y)) may be 

used as a starting point to evaluate the wave transmission coefficient at z. 
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The solution of Eq.0 can be represented through the eigenfunctions of 
Hamiltonian Eq.0 

<f>( x ,y,z) = 'Y J c n e~ lEnZ ^ n {x,y) (5) 

n 

and 

H4> n (x,y) = E n cf> n (x,y). (6) 

Finally, substitution of Eq. (J5J) into Eq.0 yields the solution of Maxwell 
equation 

V, z) = e lk ° z ^2 c n e~ lEnZ (j) n {x, y ). (7) 

n 

It follows from Eq. 0 that the local transmission amplitude of a central 
diffracted wave can be defined as 


t(x,y)= ^2 c n.e 

E n <ko 


—iE n L 




( 8 ) 


where L is the system size in the z direction. 

Before entering into a more detailed analysis of the local transmission am¬ 
plitude, let us note that if we ignore the losses and take into account that the 
metal dielectric constant in the optical region is a real large negative number, 
then: (i) the potential energy term V{x,y) = 4^ (1 — e(x,y)) in the Hamilto¬ 
nian Eq.0 is positive everywhere, (ii) correspondingly, all E n are also real and 
non-negative E n > 0 and (ii) exploiting \d 2 (j)/dz 2 \ « 2k$\d(j)/dz\ leads to the 
condition E n « 2fco- 

The central diffracted wave transmission coefficient that is measured in the 
experiment can be estimated by using the following expression 

1 


T = 


dxdy 


t(x,y) 


(9) 


( 10 ) 


where S is the area of the system. Substituting Eq.0 into Eq.0, one has 

r = 4 E M 2 - 

E n <.ko 

In order to find the coefficients c n let us consider the Eq.0 for z = 0 

$(a;, y, z = 0) = ^ y). (11) 
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Next, we assume that the impinging to the system wave has an amplitude 1 (the 
region z < 0). From the continuity at z = 0, one has $(£, y, z = 0) = 1 + r(x, y), 
where r(x, y) is a local reflection coefficient which for the metal without holes 
is approximately — 1. Clearly, the existence of the holes will change the value 
of r. However, this change will not affect the further calculations, and for this 
reason in our further calculations, for the reflection coefficient we assumed some 
average value r close zero provided that metallic fraction is small. Within this 
approach, multiplying both sides of Eq. ED by 4>n(x, y) and integrating over the 
surface, one has 



( 12 ) 


Substituting Eq. ED into Eq. ED, we arrive at the final result for the transmis¬ 
sion coefficient 



(13) 


where p = (x, y) is a two dimensional vector on the xy plane. 

This is our main general result. In the following subsections we analyze its 
limits for different models. 

To close this section let us note that if the dielectric permittivity e(x, y) is a 
periodical function, then the spectrum of Hamiltonian Eq. Jp consists of allowed 
and forbidden energy bands. As for the transmission coefficient, it depends on 
the position of incident wavenumber in the transverse energy spectrum. 

3. Bloch states 

As it follows from Eq. ED the transmission coefficient equals zero provided 
that fc 0 < Eb, where Eb is the bottom value of first energy band. When fco 
lies in the zone, using Bloch states transmission coefficient can be rewritten in 
terms of a quasi-momentum q as 


T = |l+r| 2 /, 


(14) 


where 



(15) 
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Here integration over quasi-momentum q is carried out over the first Brilloin 
zone —7r /a < q x < tt/ a, —n/b < q y < ir/b and a, b are periods of e(x,y) in the 
x and y directions, respectively. According the Bloch theorem the eigenstate cj)^ 
in a periodical potential can be represented in the form 

Mp) = ( 16 ) 

where u^(p) is a periodical function satisfying the equation 

-^W+V ) 2 + H(p) 

For simplicity and as an illustration of the approach, we will carry out further 
consideration in one-dimensional case. We hope that the results, obtained in 
this particular case, will enable us to understand, at least in the qualitative 
level, the transmission of an electromagnetic waves in three-dimensional system 
in the presence of a transverse waveguide modes tunneling. 


u ? (p) = E(q)u ? (p). (17) 


Kronig-Penney model. Suppose that slits are periodically placed in the x axis, 
which is transverse to the direction of propagation. A cross section of the 
potential in the x direction can be presented as an array of square potential 
wells. A metal part will serve as a barrier and characterized by width b and 
period a. A width of a slit is a — b, correspondingly. The metallic dielectric 
constant described by the Drude model e m = 1 — w 2 /w 2 and the height of a 
barrier is defined as V m = fc 2 /2fc 0 (k p = w p /c and uj p is the plasma frequency of 
a metal). The vacuum part dielectric describes with a e = 1 and with a potential 
energy V = 0. For a metal in optical region usually V m > k 0 . Because only 
the energies E n < ko give contribution to the transmission coefficient T we will 
consider the case E < V m when finding the spectrum of Hamiltonian Eq. d4j. 
The quantum-mechanical problem Eq.(U) is reduced to the well-known Kronig- 
Penney model 
case 


12]. Let us write /, defined by Eq. (1151) . for one-dimensional 
^ Jdxdx'rj>* q {x)ip q (x'). (18) 




E(q)<k 0 
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Using Bloch theorem (j) q (x) = e lqx u q (x ), where rt g (x)-periodical function one 
obtains 


1 r da f"“ r ma 

I= — / — / dxe~ iqx u* q {x) / dxe iqx u q (x). (19) 

-L JE(q)<ko nrn J(n—l)a J (m—l)a 

Changing the variables one finds 

/ = —— f V e - I,jan X V e lqam f dxe~ iqx u*Jx) [ dxe iqx uJx), 

L * JE( q )<k 0 2tt ^ J 0 -/o 

( 20 ) 

where wave function in the unit cell and in the different regions, is found from 

Eq.tHI 


u q i(x) = (A cos /3x + B sin f)x)e lqx , 
u q 2 {x) = (Acoshax + Dsmhax)e~ lqx , 


0 < x < a — b 
a — b < x < a 


( 21 ) 


with /3 = y/2k 0 E and a = 
Substituting Y^ n e~ mqa 


y/ 2fco(E m - E). 

= 27 rS(qa) into Ea. (l20l) . one obtains 


I = 


[ a 

/ u(x)dx 


( 22 ) 


« Ido 

where u(x) = u q= o(x) is determined by Eq. (ED- The constants B, C, D can 
be expressed by A using boundary conditions. A itself can be found from the 
normalization condition \u(x)\ 2 dx = 1/N, where N = L x /a is the number 
of unit cells. Using continuity at x = 0, one gets C = A. From the continuity 
of du/dx at x = 0, one has B/3 = Da. Finally from the continuity of u{x) and 
du/dx at x = a — b and periodicity du\{x)/dx\ x - a -b = dun/dx\ x --b , one has 


A [cos j3(a — b) — cosh ab\ + B 


sin (3 (a — b) + — sinh ab 


= 0 


A [a sinh oh — /Isin/^a — b)\ + B [/3cos/3(a — b) — /3 cosh ah] = 0. (23) 


Equalizing the determinant of 2 x 2 homogeneous equation to 0 one gets the 
dispersion relation for q = 0 

a 2 — /3 2 

1 = —-—-— sinh ab sin j3(a — b) + cosh ab cos /3(a — b). (24) 

2 ap 
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Figure 2: Plot of the function I(b/a) on a metal fraction for two values of k p (silver and 
potassium). The wave number of an incident photon ko = 10" 2 nm _1 . 


Expressing now coefficient B through A according to Eq. (El) and taking 
elementary integrals we finally arrive at the expression for I: 


1 = 


a tan ^ + /3tanh(a&/2) coshaa 


-i 2 


a/3a 


2a tan ^ a 2 ^ + 2/3 tanh ^ cosh 2 aa + aufi + a/3 tan 2 ^ ^a — b— (3 2 ba 2 


Note that above expression has been derived using the first relation of Eq. (Ell . 
However, using dispersion relation (El) , one can show easily that the second 
relations leads the same result. In figure 2 we present the function I versus b/a 
for k p = 4.6 x 10 _2 nm^ 1 (silver Ag) and k p = 1.84 x 10 _2 nm _1 (potassium 
K) (in both cases the wave number of an incident photon ko = 10 ~ 2 nm~ 1 ). 
We have checked that presented curves are unaffected by the change of ko in 
the visible region 0.62 x 10 -2 < ko < 1-57 x 10~ 2 nm _1 . This means that for 
a given period a and k p . all the curves can be scaled into a single curve. The 
physical reason of the ko independent of the results is in the structure of the 
dispersion relation (l24l) . The latter, in the frame of the adopted approach, can 
be described by a mean of the unique combination koE. 

We now turn to the numerical calculations of the I dependence on b/a for 
some typical values of a = 100 — 500nm and b = 10 — 500nm in visible region 
0.62xl0 -2 < k 0 < 1.57x 10 _2 nm _1 . The results are plotted in Fig.2. As follows 
from Fig. 2, even a small amount of metal is enough to essentially reduce I from 


.(25) 
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its maximal value 1 to almost 0 around an origin. This jump happens on a short 
scale (b/a ~ 0.01), where three curves reached their minima. The minimum can 
be explained by diffraction, which spreads and reduces the initial light intensity 
randomly across the entire system. For silver, (circle and diamond symbols) 
I increases very slowly with increasing a fraction of metal at the beginning 
and remains almost flat with further increasing of b/a. However, for potassium 
with relatively small metal barrier height and with large tunneling rate across 
a barrier, the central diffracted wave transmission coefficient (square symbol) 
increases with further increasing of a metal fraction (see Fig. 2). This leads to 
some focusing effect in a central diffraction direction. 

Note, that the transmission coefficient T, Eq. (El. includes also the multi¬ 
plier |l + rj 2 which in contrast to the I, decreases with increasing b/a. However, 
the rate of decreasing of |1 +r| 2 is much slower than the rate of increasing of I, 
and as a result the central diffracted wave transmission coefficient is increasing 
with b/a. Such an unusual behavior is caused by periodicity and Bloch states. 
It is easy to get convinced that in the case of isolated slits, where the wave func¬ 
tions overlap is negligible and where Bloch states cannot exit, the T behaves 
totally differently. 

4. Isolated slits 

For large plasmonic wavenumber k p band’s width becomes very narrow. The 
reason for it is that k p determines the barriers height and large k p suppresses 
tunneling through the barrier. In this case the transversal wave functions be¬ 
come less and less extended in space and more localized within a holes with 
negligible overlap. In this limiting case one can use the infinite potential well 
approximation to evaluate the transmission coefficient m- Writing the wave 
functions in the form 




( 26 ) 
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and substituting Eq. tE^l) into Eq. (fl3l) . we find 

m _ 8|l + r| 2 a-b 

is o 

n 


( 27 ) 


We arrived at the above expression summing over all the independent slit contri¬ 
bution and restricted ourselves by terms n = 1 while calculating sum in Eq. (1131) . 
The contributions of terms with n > 1 become irrelevant because fco < E n and 
therefore only the first band give contribution to the transmission. 

Comparing Eo. (|27| with the maximal value, one has Ti S /T max ~ As 
expected, for isolated holes the transmission coefficient becomes size dependent, 
that is Ti s /T max is proportional to the fraction of the vacuum part in the sys¬ 
tem. In this sense, the isolated holes system reveals a usual dependence of 
transmission coefficient on metal fraction. We expect to find the same ratio be 
valid also in the case of two dimensional hole array. 


5. Disordered hole arrays 


In this case it is convenient to represent the transmission coefficient, Eq. (fl3l) . 
in the form 

h -L r \ 2 r k ° 

T = ^ 1 / dE<J2HE-E n ) |<M0)| 2 > (28) 


where < ... > means averaging over random positions of holes and (f>n(o) is 
the Fourier transform of ifinif) satisfying Schrbdinger equation with random 
potential 




4>n{r) = E n (f) n (f) 


(29) 


where V (r) is assumed to be Gaussian distributed random function with a 
correlator B 


<{V(r)-V)(V(r')-V)>=B{\r-f'\) (30) 

where V = g f drV(r) = k 0 ( 1 — e m )(l — f v )/2 and f v is the fraction of the 
vacuum part (in ID periodic system, discussed in the previous subsection 1 — 
f v = b/a). 
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We now turn to the calculation of the transmission coefficient, Ea. (1551) . In 
order to carry out averaging over the randomness, it is convenient the latter 
quantity to express through the average Greens function 

fko 


T = |l + rf / — < -ImG E (q = 0) > 

Jo ^ 


(31) 


with Ge = [E — H + iS\ 1 . 

The averaged Green’s function can be represented in the form 
< G E (q) >= 1 


Q 


(32) 


E + V- £ 

where £ = YJ , n >2 the self-energy constituting contributions of irreducible 
parts of different order. In further we will restrict ourselves by the first term in 
the sum 


E 2 ($) = 


dk 

(2tt) ; 


B(\q - k\)G 0 (k) 


(33) 


where Go(q) = [E — q 2 /2k 0 + z£] -1 is the bare Green’s function. 

The explicit form of correlation function B(q ) is needed to obtain a closed 
analytical expression for T, Eq. (ED- Particularly, in the limit of very small 
size h —> 0 holes, B(q) can be substituted by B(q = 0) = B 0 ~ /„fcg/i 2 (l — e m ) 2 . 
Evaluating the integrals Eqs. (ED and (l33l) we find 


T d = 

When obtaining Eq. 
tions in the limit Bq 
we finally obtain 


|1 + rT 


arctan 


2(fcp + V ) 

k 0 B 0 


— arctan - 


2V 


koBo 


(34) 


we neglect ReE relative to V. Expanding arctan func- 
0, for the transmission coefficient in the disordered case, 


T d = 


|1 + r | ; 


k^Bo 


(35) 


2tt V(k 0 + V) 

Comparison of the Eas. (l271) and (1551) shows important differences between two 
cases, in spite of and due to the presence of the formal factor |1 +r(w)| 2 (see. 
Ref. [15(). In disordered case, Eq. (ED, in contrary to periodical case, r(w) has 
no peculiarities and is a smooth function of w. This means that the randomness 


destroy the resonant spectral shape and lead to its broadening [15|, 18(, assuming 
that the second multiplier in Eq. (1351) is smoothly varying function of the 
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frequency. This is true for almost all accepted metallic models and can be seen 
using the explicit form of the dielectric constant 1 — e m = u> 2 /uj 2 , B 0 and V. 

By comparing with isolated case contribution, Ea. (|27l) and assuming that 
V ko and B 0 —» 0, one finds 

jr ~ tfh 2 (36) 

J- is 

It follows from Eq. (1551) that disordered case the transmission coefficient Td much 
smaller than Ti S provided that koh <^. 1. Two coefficients become of the same 
order when k 0 h ~ 1. 

6. Conclusion and Discussion 

We have discussed the problem of light transport through a perforated 
metal, taking into account the transverse waveguide modes tunneling in two- 
dimensional plane perpendicular to traveling direction of light. Periodic, isolated 
and disordered holes systems are analyzed in details. Analytical expressions are 
derived for all different regimes. The one dimensional case of periodic holes, 
due to its simplicity, is investigated in detail. In the dilute metal regime, when 
metal fraction is small, our numerical study of the transmission coefficient of 
central diffracted wave indicates the existence of a minimum which is completely 
independent of an incident wavelength. The transmission coefficient of central 
diffracted wave increases when metal fraction of the system is increasing. The 
main contribution to the transmission coefficient is connected with extended 
states that close to the center q = 0 of the Brillouin zone. This means that 
in order to observe the above mentioned peculiarities in perforated systems, it 
is enough that the system exhibits long or quasi-long-range structural order in 
xoy plain (see also Ref.Q). 

In our discussion we take into account the influence of a transverse tunneling 
between different holes on the transmission coefficient T. As a result, T does 
not depend on the system thickness in z direction(we ignore the imaginary part 
of e). Obviously, exponential decaying of T with thickness will arise on if one 
takes into account the losses. 
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In our consideration we substitute a local reflection coefficient by an average 
value. This assumption seems more relevant in the random hole arrays case. 
However, even in the periodical case one can imagine r(w) as a quantity that 
correctly takes into account periodical hole arrays similar to periodical gratings 
[l7|. Note that all plasmonic effects are included into r(w). Particularly on 
impinging of p-polarized light plasmon is generated on the perforated surface. 
Reflection coefficient close to the plasmonic resonance becomes minimal 113 


leading to the maximal value of transmission coefficient. Plasmonic resonance 
takes place when plasmon wavenumber coincides with one of the photonic crystal 
reciprocal lattice periods,see for example, 0. Note that r(w) depends on the 
geometry of the perforated surface and can detect resonance effects associated 


with the geometry of holes 


H- 
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